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| J . . Jt  has  1)6611  known  for  nearly  a century  [1]  that  a rotating  column  of  fluid  can  ex- 

I hlbit  a class  of  neutral  modes  called  inertial  oscillations.  It  is  plausible  to  expect  concen- 

trated atmospheric  vortices  such  as  tornadoes,  water  spouts,  and  dust  devils  to  exhibit 
inertial  oscillations  under  proper  conditions  of  forcing.  As  convection  becomes  established 
in  a large  cumulonimbus,  surrounding  low-level  air  is  sucked  into  the  base.  This  air  will 
approximately  conserve  potential  vorticity  as  it  converges,  thus  causing  the  initiation  of  a 
concentrated  vortex.  As  the  convection  accelerates,  the  rotating  updraft  may  overshoot 
its  equilibrium  level  more  or  less  throughout  the  cell,  causing  a readjustment  toward 
equilibrium  in  the  form  of  a general  downward  motion  of  air  in  the  cell.  This  situation 
would  adjust  to  equilibrium  through  one  or  more  damped  inertial  oscillations.  The  up- 
draft portions  of  these  oscillations  have  an  accompanying  strong  horizontal  convergence 
in  the  lower  levels  over  and  above  that  associated  with  the  initial  inflow,  which  in  turn 
can  produce  the  intense  vortex  observed  as  a tornado,  water  spout,  or  dust  devil.  As  the 
updraft  subsides  and  the  downdraft  sets  in,  a low-level  outflow  results  which  diverges 
angular  momentum,  thus  causing  the  vortex  to  weaken.  This  mechanism  would  require 
that  the  life  span  of  such  a vortex  should  be  closely  related  to  the  inertial  period  for  it. 

In  fact  the  inertial  period  should  be  a lower  bound  on  the  life  cycle.  In  many  cases  tor- 
nadoes are  observed  to  occur  multiply  from  one  long  lived  parent  cell  [2-4] . In  this  case 
the  interval  of  occurrence  of  the  tornadoes  should  be  strongly  controlled  by  the  inertial 
period  for  the  rotating  system,  if  indeed  such  a situation  is  represented  by  a series  of 
inertial  oscillations. 


The  suggestion  of  a near-regular  vortex  oscillation  was  spurred  on  by  some  simple 
laboratory  experiments  in  which  a vortex  was  generated  in  a rotating  tank  of  water.  The 
vortex  was  always  oscillatory,  and  during  the  updraft  portion  of  the  oscillation  a strong 
concentration  of  angular  momentum  occurred  at  the  base,  as  one  would  observe  during 
tornado  formation.  The  oscillation  period  was  adequately  given  by  the  inertial  oscillation 
period  when  the  core  was  treated  as  a cylinder  in  solid  rotation.  This  prompted  the  calcu- 
lation of  inertial  periods  for  atmospheric  vortices  (combined  Rankine  vortices),  with  the 
result  that  the  periods  were  generally  a good  lower  bound  on  the  life  span  of  these 
vortices.  In  the  case  of  multiple  tornadoes  their  occurrence  interval  was  somewhat  longer 
than  their  inertial  period.  Thus  it  is  felt  that  inertial  oscillations  of  atmospheric  vortices 
occur  and  exert  a strong  dynamical  influence  on  the  vortices  when  the  external  forcing  is 
suitable. 
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Another  major  impetus  for  this  report  is  that  several  investigators  who  have  simu- 
lated laboratory  vortices  [5-7]  have  generally  reported  vortices  which  are  steady  or  have 
rapid  transitory  perturbations.  In  this  latter  case  no  attempt  has  been  made  to  identify 
the  oscillatory  phenomena.  In  the  present  experiments  the  oscillations  could  not  be 
eliminated,  no  niatter  what  precautions  were  taken  to  do  so,  which  indicates  they  are  an 
important  feature  of  a vortex  which  is  not  directly  coupled  to  the  force  which  generates 
the  updraft. 


LABORATORY  EXPERIMENTS 

The  vortices  were  simulated  in  a large  Pyrex  cylinder  mounted  on  a variable-speed 
turntable  (Fig.  1).  Water  was  circulated  through  a pump  and  flowmeter  and  then  upward 
from  a nozzle  in  the  center  of  the  cylinder  near  the  top.  This  water  flowed  directly  out 
to  the  wall  of  the  cylinder  near  the  top,  where  it  was  removed  with  a manifold.  The  flow 
of  water  across  the  top  induced  a secondary  flow  across  the  top,  down  the  sides,  in 
across  the  bottom  and  finally  up  the  vortex  along  the  axis  of  the  cylinder.  The  flow 
above  the  nozzle  was  turbulent  for  all  the  flow  rates  used.  To  make  the  flow  visible,  ink 
was  injected  just  above  the  base  of  the  cylinder  and  near  the  center. 


Fig.  1— Variable-speed  rotating  cylinder,  the  pri- 
mary water-flow  circuit,  and  the  induced 
secondary  flow  which  produces  the  central 
vortex 
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The  quantities  varied  during  the  measurements  were  the  flow  rate  through  the  nozzle 
Q and  the  rotation  rate  of  the  turntable  St.  For  all  combinations  of  Q (5  to  25  cm3  s'1) 
and  St  (0.25  to  3.0  rad  s'1 ),  a concentrated  vortex  was  observed  in  the  center  of  the 
cylinder.  The  vortex  appeared  to  have  a double  structure— a concentrated  core  which  had 
a large  angular  velocity  cu„  and  a wall  about  the  core  which  rotated  at  a rate  larger  than 
SI  but  less  than  ov  This  structure  is  clearly  evident  in  Fig.  2.  The  core  is  readily  discern- 
ible at  the  top  and  near  the  bottom  of  the  vortex.  Both  the  core  and  the  wall  oscillate 
vertically,  but  the  net  flow  is  upward  in  both. 


In 

Fig.  2— Central  core  and  outer  wall  of  the  vortex.  The 
ink  is  injected  from  hypodermic  tubing.  Inflow  to  the 
vortex  at  the  base  has  just  been  completed. 


Typically  the  vortex  was  3 cm  in  diameter,  with  some  minor  variation  when  Q and 
SI  were  varied.  For  the  range  of  variables  used  in  the  measuerments,  ojv  varied  from 
1.5J2  at  low  Q up  to  17.5J2  at  the  largest  Q.  There  was  only  a weak  dependence  of 
(jjy /St  on  St. 

In  the  vortex  oscillation  the  center  of  the  vortex  moved  up  or  down  while  the  fluid 
immediately  adjacent  to  the  core  moved  in  the  opposite  sense  relative  to  the  net  upward 
flow  in  the  vortex.  A pure  standing  mode  in  the  vertical  was  not  observed.  Generally 
there  appeared  to  be  a mixture  of  the  first  several  modes  occurring.  A sequence  of  photos 
from  a motion-picture  record  at  St  = 1.5  rad  s’1  and  0 = 15  cm3  s*1  is  given  in  Fig.  3; 
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the  frames  are  1/4  revolution  of  the  apparatus  apart.  The  ink  (the  dark  conical  area 
which  develops  in  the  center)  is  injected  about  2 mm  above  the  base  and  about  2 vortex- 
core  radii  from  the  center  of  the  cylinder.  At  frame  4 the  core  of  the  vortex  has  reached 
its  maximum  upward  excursion  and  maximum  coy.  It  then  collapses  downward  through 
frame  6,  and  decreases  markedly  with  the  collapse.  The  core  then  rises  until  frame  11, 
collapses  slightly  and  grows  again  until  frame  17,  when  it  collapses  again.  Another  col- 
lapse occurs  after  frame  23.  The  double  structure  of  the  vortex  is  evident  in  frames  26  to 


The  periods  of  these  oscillations  were  measured  over  a fairly  wide  range  of  condi- 
tions. The  periods  are  not  highly  regular,  presumably  due  to  the  mixture  of  modes  pre- 
sent, so  averages  were  used  in  most  cases.  In  Fig.  4 the  angular  frequency  of  the 
oscillation  0iOsc>  measured  visually  with  a stopwatch,  is  plotted  against  the  average  rota- 
taion  rate  of  the  vortex  core,  also  measured  visually  with  a stopwatch  (and  is  compared 
with  lines  calculated  using  equations  not  yet  discussed).  The  vertical  bars  show  the  range 
of  frequencies  observed. 


Fig.  4— Comparison  of  measured  and  calculated  periods.  The  solid 
lines  are  from  Eq.  (1)  and  the  dashed  lines  are  from  Eq.  (3)  for 
00  = 2.405.  The  range  of  measured  points  is  due  to  the  n irregularity 
of  the  oscillations. 
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The  frequency  wOJC  of  a standing  inertial  mode  having  n - 1 nodal  surfaces  in  the 
vertical  (excluding  the  ends)  in  a cylinder  of  height  d and  radius  ro  in  solid  rotation  is, 
according  to  Fultz  (8], 
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the  latter  being  so  when  X,d  > n X r0.  Here  X,  is  the  ith  root  of  J\  (X,)  = 0 where  J\  is 
the  ordinary  Bessel  function  of  the  first  order.  For  the  present  geometry  the  curves  for 
n = 1,  2,  and  3 (X,  = 3.833)  are  drawn  in  Fig.  4 as  the  solid  lines.  The  value  of  r0  used 
is  1.5  cm,  which  was  obtained  from  the  photographs  of  Fig.  3,  and  d was  taken  as  39 
cm,  and  effective  depth.  Typically  the  dominant  n observed  was  2 or  3.  The  measured 
frequencies,  although  somewhat  irregular,  do  fall  into  the  inertial  range;  hence  it  is  con- 
cluded the  observed  oscillations  are  mixtures  of  several  modes  of  inertial  oscillations. 

The  source  of  the  oscillations  is  obscure;  however,  since  inertial  oscillations  are 
neutral,  they  will  occur  at  or  near  their  natural  frequency  when  any  forcing  function  is 
present  which  has  components  near  the  inertial  frequencies.  The  small  variations  in  the 
rotational  drive,  variations  in  the  pump  flow  rate,  and  instabilities  in  the  nozzle-manifold 
flow  all  contribute  to  the  forcing.  The  important  facts  are:  the  fluid  oscillations  had  the 
correct  structure  and  roughly  the  proper  periods  so  they  could  only  be  inertial  oscilla- 
tions. 


THEORY 


To  approximate  the  flow  field  of  an  atmospheric  vortex  more  accurately  than  the 
theory  which  led  to  Eq.  (1),  i.e.  homogeneous  fluid  in  solid  rotation  in  a finite  cylindri- 
cal container,  a basic  flow  field  which  is  a combined  Rankine  vortex  in  a stratified  atmos- 
phere will  be  used.  There  is  ample  evidence  that  such  a flow  field  exists  in  tornadoes  [9] 
and  dust  devils  [10].  It  will  further  be  assumed  this  basic  flow  field  is  horizontal  and 
height  independent.  This  latter  assumption  is  quite  crude,  but  the  intent  here  is  to  obtain 
approximate  oscillation  periods.  The  small-amplitude  theory  will  be  used;  it  is  quite 
reliable  for  estimating  finite-amplitude  periods  although  generally  inaccurate  in  specifying 
any  more  than  the  gross  features  of  the  flow. 


The  tangential  velocity  in  the  combined  Rankine  vortex  is 

v = rnc,  0 <r  <rc, 

= ncr*/r,  rc<r<» 


(2a) 

(2b) 


where  ftc  is  the  core  rotation  rate  and  rc  is  the  core  radius.  The  oscillation  period  is  then 
(the  details  are  left  to  Appendix  A) 


(3) 


where 


and 


VT  rcf2c,  J0(/3q)  - 0, 


(4) 
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in  which  H is  the  height  of  the  tornado  vortex  and  of1  is  the  equivalent  depth  of  the 
atmosphere  when  the  density  is  given  by  p = p0e-<“  (a'1  = 9.5  X 103  m). 

It  is  assumed  for  (3)  that  y < 1,  % < y2r$t  aV2/g  ^ j For  yr  = 100  m rl 
this  term  is  0.1,  and  this  value  is  approached  only  for  large  or  intense  tornadoes.  This 
result  is  similar  to  (1),  except  here  the  roots  of  J0  and  the  factor  7'1  enter. 

The  vertical  velocity  field,  obtained  from  (.44)  and  (.49),  is 

W = ~)eaz  12  cos  at,  0<r<rc,  ( 


~ yAK1(yrc)  Ko(yr)eazlz  sin  — cos  at,  rc<r<  oo  . (5b) 


COMPARISON  WITH  OBSERVATIONS 

The  explicit  value  of  /3q  has  not  been  chosen  as  yet.  The  experiments  showed  two 
distinct  regions  of  vertical  motion  in  the  vortex,  but  it  was  not  clear  if  they  were  some- 
times opposite  in  direction  because  of  the  net  vertical  flow  superimposed  on  the  vortex. 

A similar  double  structure  was  observed  by  Sinclair  [10]  in  dust  devils.  Sinclair’s  data 
shows  a definite  reversal  in  w.  From  (5)  this  occurs  when  fig  is  chosen  as  the  second  zero 
of  Jo(Po  - 5.507).  Some  tornadoes,  such  as  the  Dallas  tornado,  show  a strong  tendency 
toward  a first-mode  structure.  Then  p0  = 2.405. 

According  to  (3)  the  period  of  the  oscillation  requires  a knowledge  of  VT  and  H 
(since  y is  a function  of  H).  We  are  assuming  r represents  either  the  duration  of  the 
phenomena  or  the  occurrence  interval  in  the  case  of  successive  multiple  tornadoes  being 
spawned  from  the  same  parent  cloud.  Until  recently  data  on  Vy  and  H were  almost  im- 
possible to  obtain.  Burgess,  Lemon,  and  Brown  [11]  have  recently  obtained  doppler  radar 
data  from  a tornado  which  touched  down  at  Union  City,  Oklahoma,  and  produced  intense 
damage.  They  found  a region  of  intense  shear  in  the  cloud  with  an  average  maximum 
tangential  velocity  of  27  m s'1  at  0.2-km  and  3.5-km  elevations.  They  state  that  this 
intense  shear  zone  (presumably  the  tornado  vortex)  extended  to  10  km  and  retained 
vertical  continuity  through  this  range.  The  cloudtops  on  radar  were  at  a height  of  15  km. 
The  strong  shear  region  observed  near  the  ground  level  definitely  corresponded  with  the 
observed  tornado  position.  This  confirms  that  a tornado  is  a continuous  vortex  extending 
vertically  through  most  of  the  cloud  mass. 

The  tangential  velocities  observed  in  the  cloud  are  much  less  than  those  which  pro- 
duced “intense  damage”  at  ground  level  (which  are  denoted  emax,  d)-  Based  on  Ref.  4, 
considerable  damage  occurs  for  winds  of  50  to  70  m s'1  (F2  winds),  and  severe  damage 
occurs  for  winds  of  70  to  90  m s"1  (F3  winds).  The  Union  City  tornado  had  a transla- 
tional speed  v of  8 knots,  so  the  “damage”  tangential  velocity  (emax  d - v')  was  62  to 
82  m s 1 , or  1.5  to  3 times  the  tangential  velocity  observed  in  the  cloud. 
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Also,  if  (3)  is  used  to  estimate  periods,  Vj  represents  an  average  tangential  velocity 
whereas  the  velocity  observed  by  Burgess,  Lemon,  and  Brown  [11]  would  be  a maximum  ’ 
value  Umax , because  this  is  when  the  vortex  reaches  sufficient  intensity  to  manifest  itself 
If  we  let  a represent  the  amplitude  of  the  inertial  oscillation  as  a fraction  of  the  basic 
state,  then 


1 + “ ° %a */VT- 


(6) 


Also,  if  s represents  the  ratio  of  the  damage  tangential  velocity  (umax  A - v ) to 
then  (3)  becomes  ’ max’ 


r =- 


7r0o(l  + a) 


(7) 


or 


r =- 


7T(30s(1  + a) 


max,  d ~ ^ 


(8) 


The  Union  City  tornado  suggests  s =»  2.  The  amplitude  a is  uncertain,  but  a range  of  0.2 
to  0.5  is  reasonable.  The  basic  state  velocity  VT  for  our  theory  should  be  the  tangential 
velocity  of  the  parent  rotating  updraft.  This  is  obtained  here  through  the  factors  a and  s. 


Data  for  several  documented  tornado  occurrences  (Dallas  [4,  9] , Fargo  and  Palm 
Sunday  [9],  and  Union  City  [11])  along  with  dust-devil  data  [10]  and  waterspout  data 
[Golden,  private  communication]  are  used  for  comparison.  For  the  Fargo  and  Palm  Sun- 
day occurrences  umaXi  A is  available.  For  Dallas  the  largest  velocities  Koecker  [9]  gives 
are  from  30  m to  about  250  m above  the  ground  and  are  75  m s'1.  The  damage  esti- 
mates from  Dallas  give  velocities  ranging  from  45  to  80  m s'1 . The  cloudtop  height  is 
used  for  H.  For  Dallas  the  tropopause  height  was  used,  and  for  the  Palm  Sunday  torna- 
does 15  km  was  arbitrarily  chosen.  The  occurrence  interval  for  the  tornadoes  is  taken 
from  Fujita  [4],  For  the  Fargo  and  Palm  Sunday  tornadoes  Eq.  (8)  was  used,  and  for  the 
other  cases  Eq.  (7)  was  used. 


The  values  of  r for  a = 0.2  and  0.5  are  given  in  Table  1.  These  calculated  periods 
are  also  plotted  against  the  observed  periods  in  Fig.  5,  differentiated  by  mode  and 
whether  a duration  or  interval  of  occurrence  was  given.  The  inertial  period  is  more  closely 
related  to  the  interval  than  to  the  duration;  however  the  duration  should  be  a significant 
fraction  of  the  period,  depending  on  the  swirl  intensity  of  the  parent  storm.  In  any  case 
the  duration  should  be  bounded  above  by  the  inertial  period. 


CONCLUDING  REMARKS 


This  calculation  of  the  inertial  period  neglects  many  factors,  the  most  important  of 
which  are  nonlinear  effects  and  interactions  of  the  vortex  with  the  ground.  It  is  intended 
to  show,  within  a factor  of  2 over  a period  range  of  two  decades,  that  inertial  periods 
agree  with  natural  vortex  lifetimes  or  occurrence  intervals  and  that  the  presence  of  an 
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Fig.  5— Comparison  of  calculated  and  observed  periods  for 
dust  devils,  waterspouts,  and  tornadoes.  Data  for  the  duration 
of  the  vortex  and  for  the  interval  between  occurrences  are 
given.  Both  modes  are  included. 


inertial  oscillation  on  a vortex  should  be  expected  when  variable  natural  forcing  is  pre- 
sent. The  inertial  oscillation  provides  a mechanism  whereby  angular  momentum  can  be 
concentrated  into  an  intense  vortex.  The  use  of  a linear  theory  here  is  an  approximate 
treatment  of  the  problem.  Numerical  or  other  nonlinear  time-dependent  calculations 
should  provide  more  accurate  estimates  of  intervals  or  durations,  if  inertia  oscillations  are 
the  fundamental  mechanism  controlling  atmospheric  vortex  formation. 
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Appendix  A 

DERIVATION  OF  EQ.  (3) 


Consider  a cylindrical  coordinate  system  in  which  r,  0,  and  z are  respectively  the 
radial,  tangential,  and  vertical  coordinates.  The  oscillation  velocity  components  in  the  r, 

6,  and  2 directions  respectively  are  u,  v,  and  w.  Then  the  appropriate  linearized  equations 
of  motion  for  axisymmetric  motion,  with  a basic  axisymmetric  velocity  V = V(r),  are 


( bw\  bp 

p\-bi)  " -17’ 


(Ala) 


(Alb) 


(Ale) 


bu  u bw 

5-  +—  + = 0 

or  r bz 


(Aid) 


(Ale) 


n dp 

0 * - 37  - *>•  <aih 

where  P is  the  pressure  distribution  for  the  basic  state,  p is  the  perturbation  pressure,  and 
g is  gravity.  Since  the  atmosphere  is  stratified,  let 

P = Po(r)e'“.  (A2) 

where  a'1  = 9.5  X 103  m.  From  (Ale)  and  ( Alf) 

±dpo  =aV>  (A3) 

Po  dr  rg 


jj 
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Introduce  a Stokes’  stream  function  tp  as 


1 dtp 

u — - 

r dz  ’ 


1 dtp 

w = 

r or 


so  (Al)  becomes  upon  elimination  of  u , v,  w and  retention  of  the  r and  z density  varia- 
tions 

H a /i a^\]  /a  \r a2  2v/v  avYiiag. 

l -*«  * Sr  t */J  fe  ' 7U72 ' 7l;  * ZKPZ  ‘ ° ' (A! 


The  solution  is  of  the  form 


tp  = i//(r)(cos  ot)e?*/2 sin  ^ , (A( 

where  o is  a constant  to  be  determined  and  H is  the  total  height  of  the  vortex.  Substitution 
of  (A6)  into  (A5)  yields 

3 2tP'  faV2  \1  dp'  [2  V (V  dV\  9 , 

* It-1);*  *b(r7)-r  = o'  <ai 


in  which 


2 7I'2  a2 

r = — + — . 

H 2 4 


The  cdmbined  Rankine  vortex  is  described  by 

V = rSlc,  0«Sr<rc  (inner  region,  I),  (A8: 

= rc&c/r’  rc<:r<:°°  (outer  region,  O).  (A81 

In  both  regions  (A7)  has  solutions  in  closed  form  when  aV2/g  = 0.  For  the  most 
intense  tornadoes  V < 100  m s"1;  hence 

— < 10'1, 

g 

and  this  term  can  be  neglected  in  comparison  to  unity.  Then  the  solutions  for  the  inner 
and  the  outer  regions  are  respectively 


